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For Gagen-Poiseuille flow, we show that exponential instability (to extremely small, axially 
symmetric disturbances represented by Galerkin's approximation) is possible only if there exists 
bi-periodic variability of the disturbances along the pipe axis when the threshold Reynolds 
number Re t iXp) depends on the ratio p of two longitudinal periods. Absolute minimum 

Re.7(j?)«448 (for /?«1.53..) is obtained that corresponds to the observed conditions of 

transition from the laminar resistance law to the turbulent one and Tollmien-Schlichting waves 
exciting in the boundary layer. 
PACS: 47.20.Ft, 47.27.Cn, 47.27.nf 

1. Fundamental and applied problem of understanding of the mechanism of turbulence 
development for the Gagen-Poiseuille (GP) l flow remains mysterious because of the paradox of 
the linear stability of this flow to extremely small magnitude disturbances for any large Reynolds 
number values Re [1-4]. Up to now, to bypass this obvious paradoxical contradiction to 
experiment, a consensus is made to allow only non-linear mechanism of the GP flow instability 
to a disturbance with a sufficiently large finite magnitude [5-8]. This assumption is usual (see [3, 
4]) based on the very special interpretation of the experiments in which the substantial increase 
of the threshold Reynolds number Re^ value (defining transition from the laminar to turbulent 

flow state) is achieved by increasing the smoothness of the streamlined pipe surface. In such an 
interpretation, only correlation between smoothness increasing and resulting decreasing of the 
initial disturbance magnitude is taken into account. Meanwhile, noted already by O. Reynolds 
[1] extremely high sensitivity of the value of Re^, to the initial disturbance does not exclude 

possibility of affecting on Re^ not only magnitude but frequency characteristics of the 

disturbances (resulting from non-ideal smoothness of the streamlined surface) as well. Actually, 
for example, in the experiment [9], it was found that under the fixed magnitude of artificially 
excited disturbances, GP flow instability develops only in the narrow disturbances frequency 
range. 

Herein, we show that possibility of the linear absolute instability of the GP flow in the general 
case also is defined by the value of the additional to Re control parameter p , which 



1 GP flow is a laminar, stationary flow of the uniform viscous fluid along the static, direct and the length unbounded 
pipe with the round the same along the pipe axis crosscut. 

2 For special GP flow modifications (in the pipe with elliptic crosscut [10], or in the rotating pipe [11], or in the case 
of the flow with the finite size particles transported by it[12]), there initially exists natural additional to Re control 
parameter, and the paradox of linear stability is absent. 
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characterizes frequency properties of the disturbances and affects on the value of Re w (p) 

independently from the magnitude of the initial disturbances. Such a parameter p introducing is 
performed below on the base of pointed by O. Reynolds [1] (and by W. Heisenberg, see [4, 6]) 
dissipative GP flow instability mechanism related with the action of molecular viscosity v in 
the proximity of the solid boundary. According to [1], the mechanism manifests itself as a 
spontaneous one-stage appearance for Re > Re ( i of the vortexes having character size / , 

«.. which does not grow any more contrary to expectations with the growth of the velocity 
magnitude [1]». As a result, in addition to the character disturbance scale related to the pipe 
diameter, 2R, due to near-boundary vortex-generating viscosity action, there exists a new 
additional scale / ,which is proportional to the size of some part (which is really activated only 

for v^O) of the pipe boundary. It, together with R , can define frequency parameters of the 
initial disturbances, e.g., the longitudinal along the pipe axis (axis z ) spatial periods. The ratio of 

/ 
the periods, p = — , as it is shown below, is a new additional to Re parameter defining the GP 
7.R 

flow instability threshold to the conditionally periodic on z extremely small vortex disturbances. 
Such representation of the disturbances structure meets in observed conditionally periodic 
Tollmien-Schlichting (TS) waves development of which (caused also by near-boundary action of 
the molecular viscosity) precedes explosive emerging of the turbulence in the near-boundary 
layer [17-19]. Moreover, already in [2, 20], it was noted that usually considered in the linear 
stability theory strictly periodic on z structures of the disturbances field obviously do not match 
those observed in the experiments. We show further, that the assumption of the strict 
longitudinal (by z ) periodicity of the disturbances may result in the paradoxical conclusion on 
linear stability of the GP flow for any Re with Re.r ->oo. The same time, with the simplest 

Galerkin's approximation , we get the finite value of the absolute minimum Re^, ~ 448 (for 

p «1.53..), that is the observed characteristic threshold value Re for transition from the laminar 

resistance law to the turbulent one [2, 21] and for the conditions of the TS waves exciting [19]. 

Qualitatively the proposed theory conclusions agree with the experimental data for the flows in a 

pipe [22-24]; also, we have conducted comparison of our conclusions with conclusions of the 

stability theory (Schlichting and Lin's), and with experimental data on the stability of laminar 

near-boundary layer [25]. 

2. Let's consider the well known (see [4]) representation of the GP flow in the cylindrical 

r 2 ff2 dpr> 

coordinate system (z,r,<p) : V Qr = V Q(p = 0,V Qz = V rmx (1 - —^) , where Fmax = gz » fluid 

dpo 
density p — const , -^- is a constant value of the pressure p^ gradient along the axis of the 

pipe with radius R , and v is the fluid kinematic viscosity coefficient. 

In the axially symmetric case (i.e. for extremely small disturbances which are not depending on 

<p), linear instability of the GP flow may be defined only by the value of the tangential 



3 Such a mechanism is realized in the systems having disturbances with the negative energy [13-16] as , for the 
threshold generation of the vortexes (rotons) in the stream of the superfluid helium in capillary [13]. 

Galerkin approximation is used to represent disturbance dependency not only on radial but on the longitudinal 
coordinate z as well when for two synchronous radial modes having two different variability periods on z the order 
of the exponential temporal evolution is one and the same. The latter causes necessity of the use of Galerkin 
approximatiom for describing disturbance dependency on z contrary to the conventional problem stating [2, 3] in the 
linear theory of the GP flow stability where independent modes are considered - each with its own order of the 
exponential temporal evolution. Appendix B contains detailed example illustrating distinctions of our linear 
stability problem stating from the conventional one (see, e.g. [2, 3]). 
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component of the disturbance velocity Vm , which before coming to the non-linear evolution 
stage does not depend on other disturbance fields and satisfies the following equation: 

dV dV V 

-zr + V 0z (r)^=v(AV v -^), (1) 

ot oz r 

where A is the Laplace operator. It follows from (1) that exponential growth of V with time, 

not depending on the reference system choice (i.e. absolute instability of the flow [2, 4]), 
evidently is not realizable as for v = so for the case of strictly periodic or localized by z 
functions Vm in (1). It is easy to verify the latter, if to multiply (1) by Vm and integrate the both 

sides of (1) by z in the limits of the period of the function Vm (or from -oo to +oo in the case 

of the localized by z function Vm ), when the contribution of the second term in the left hand 

side of (1) is zero. 

Thus far, only for conditionally periodic by z (see [2]) functions and only for v^Oin (1) it is 
possible to expect realization of the absolute dissipative instability of the GP flow for above the 

max 
threshold Reynolds numbers Re > Re .. , where Re = . 

Let's find solution of (1) in Galerkin's approximation (see [3]), describing dependence of Vm 
on r as follows 

V v =V^A n {z,t)J x {j x /-A , (2) 

where Vm in (2) for any A n meets necessary boundary conditions on r ( V m < <x> for r — , and 
Vm =0for r = R), because J, is the first-order Bessel function, and /-, are its zeroes, i.e. 
J,(/i ) = 0for any n = l,2,.. For coefficients A n , from (1), (2), we get the following system of 
equations in the dimensionless form: 

^+£.4. - ^ +«*£/•_ ^=o. (3) 

OT OX n= \ OX 

where m = \,2,..fl,T = —^,x = -p . In (3), constant coefficients, P nm , are of the form 



R^ 



R 



P nm =-^\dyy(\-y 2 )Uk n y)Uk m y) , (4) 

where J 9 is the second-order Bessel function, and the linear on y term under the integral yields 

{\ fi = ffi 
' . For N — 1 in (3), the last 

term may be excluded by the Galilean transformation, and hence for N-l, there is no 
possibility for the absolute instability of the GP flow. In relation with this, we shall further 
consider (3) for the simplest non-trivial case, N — 2, that allows already to eliminate the GP flow 
linear stability paradox, and leads to the conclusions quantitatively agreeing with the 
experimental data (see section 4) [22, 24]. 

3. Observed in the experiment field structures do not agree with the strictly periodic along the 
pipe axis disturbances (see herein above, and [2, 20]). Moreover, in [20], it is noted that different 
radial modes (defining dependency on the radial coordinate) correspond to different variability 
periods along the pipe axis. The latter can be modeled with the use in (3) of different periods 



along the pipe axis for the modes with different index value m. Such a requirement corresponds 
to introducing for each of these modes of its own, independent from other modes, periodical 
boundary condition. Hence, it is necessary to use the Galerkin's approximation for a dependency 
of the solution of (3) on x relevant to the mentioned above boundary conditions. Let, thus far, in 
(3), for example for N — 2, magnitudes Ay and A 2 have the form of the running waves with 

different periods along the pipe axis: 

A [ = A [ Qe 1 T+i2mxx ? ^ = A ^ t+ilnftx ? (5) 

where Ar>and A^q are constants, and the additional to Re control parameter is defined as 

a l 

P = -n for all a and /? , when, in particular, p = — = -^— . In Galerkin's approximation (for the 

solution of the system (3) depending on x in the form (5) for a =£ J3 generally), we get a uniform 
system for A q and Ajq , condition of solvability of which for the finite A q and Ajq allows to 

define the value of the exponent A = A l + iX 2 depending on dimensionless parameters Re , p , 
and B (see (A.l) and (A.2) in Appendix). 

Vortex structure of the considered in (2), (5) disturbances is related to non-zero components 
(longitudinal and radial) of the vorticity field, and zero value of the integral value of the Lamb's 
vortex momentum density ,which for the main undisturbed GP flow has longitudinal component 

pV 
only equal to max . This structure has finite value of the integral helicity for the resulting flow 

which is a linear superposition of the main GP flow (having non-zero the longitudinal 
component of the velocity field and the azimuthal component of the vorticity only) and of the 
disturbance field. 

Usage according to [20] of the form (5) to get solution of the system (3) just in the Galerkin's 
approximation (instead of usually applied exact solution of the system (3) for the radial modes 
having one and the same longitudinal period) may be also interpreted as a result of the possible 
modeling of the always existing stochastic heterogeneities of the pipe internal surface. 

Linear exponential instability condition with A. > Oin (5) has the form of (A.3). For Re »1, 

it may be reduced to more simple (A.4).Then for (3 from (A.5), minimizing Re , in (A.4), it 

defines the following linear exponential instability condition for GP flow (providing a lower 
boundary for an exact value of Re , , determined from (A.3)): 

Re>Re f/i ^ 2 "-^ . (6) 

P P n I Si 

n \ S 
where S = sin Tip sin — sin n{ — h/?), B = (pP.. -P„), 

p p | S | [ [ zz 

F = (r^ 2 +r^+p 2 )A 2 +(yf 2 -yI x ){\-p 2 )b 2 +2ab{ Y 1 2 -p 2 yI x ), 



~> ctgn' " ' 

n\\- P y 



4SP u P 21 p^ctg7v(p + -) 
A -B — for P. 1 , P^^ P.., P.. from (4) , for arbitrary value of p 

2 
and for A > . 
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Only for A < 0, it is necessary to use instead of (A.4) and (6) the general condition (A.3), for 

which A > in (5). 

In (6), infinite value Re , — > <x> takes place only for those finite values of p for which S = , 



I k + \±^(k + \) 2 -4 
i.e. for p- p, =k,p = p 1 = — , and p = p , = for all integer 

ft _ ft" vft+1 2 

k 

k = l,2,... For p, being inside the intervals of changing of p between any two consecutive 

values of p , , p 1 , /? , , value of Re , in (6) is a function of p having one local minimum 

ft j_ Vk" + 1 "* 

inside each of the pointed out intervals (see Fig. l,b). So far, the absolute minimum value, 
Re^ « 442 , in (6) is reached for p « 1.53.., close to the "golden" ratio p = -^— * 1.618 .. 

(i.e. the limit of the infinite sequence of ratios of two consequent Fibonacci numbers 1, 2, 3, 5, 8, 
13, 21, ..). For the same p , from the exact condition (A.3), we get close to it absolute minimum 

value Re , « 448 (see also Table in Appendix where conclusions based on (A.3) and (6) are 

compared). 

4. Stated value Re, ~ 448 corresponds to the interval of Re e 3 00-^500 no ted in the 

experimental observation of the threshold transition from the laminar resistance law (for the flow 
in the pipe) to the turbulent one [2, 21], and of the TS waves exciting in the near-wall region of 
the boundary layer [19]. Observed in [1] and other experiments (see references in [23, 24]) 

hyper-sensitivity of the value of Re , to the initial disturbances, actually, corresponds to the 

obtained in (6) dependency of Re , on p when, for example, Re , in (6) changes 

approximately 600 times with the change of p from 0,1186 to 0,1118 only. Neighboring local 

minima of Re , in (6) may also significantly differ each from the other, when for p « 2.23 , we 

have in (6) Re , ~ 1982, and for p « 3.86, we already have Re , ~ 84634. Enlarged fragments 

of the neutral stability curve corresponding to the condition (6) (see Fig. 1 , b), are represented in 
Fig. 1, a, as dependence of the value of l/2p on Re. They are plotted on the taken from the paper 
[25] figure (see Fig. 12 in [25]) representing theoretical (Lin's and Schlichting's) neutral stability 
curves and corresponding to them experimental data related to the determining of the instability 
threshold in the boundary layer. 

Considered in (2), (5) conditionally periodic on z structures of the initial disturbance field 
V qualitatively agree with observed in [23] wave-like changes of V . The latter are especially 

clearly seen in the proximity of the turbulence decay threshold for Re ~ 1750 when the ratio of 

8 13 

visible in Fig. 5d) from [23] character periods with p « — and p « — is close to the pointed out 

5 8 

above "golden" p . In Fig. 2a), taken from [24] (see Fig. 4 b [24]), it is giving the dependency 

on Re for the constant velocity of the rear front of the observed in the flow in the pipe [24] 
turbulent spot together with the experimental data [22] related to the given in Fig. 2a) interval of 
the values of Re (dark triangles for the rear front of the turbulent spot, and white triangles for the 
leading front). Also, results following from the non-linear theory [8] and for the phase velocity 



X n V V 

y = — z max ( scaled by the mean flow GP ve i oc ity y = _max_ ^ defined in ( A 6 ) from the 

P 2^ Re m 2 

representation of /Lin (A.2) for the neutral curve (i.e. when /t = 0) are presented there. 

Hence, a theory based on the simplest Galerkin's approximation with TV = 2 in (3) already 
yields agreeing with experiments [22, 24]. The same time, its worth to mention that linear 
exponential instability in (6) turns out to be possible also for rational p not equal to p = p k or 

p 1 , but for irrational p = p . , k = 1,2,.., vise versa, Re , — > qo in (6). It means that these 

~k 
conclusions based on the Galerkin's approximation with TV = 2 may in this respect differ from 
the made above (in the section 2) general conclusions on the base of the original equation (1). 
Future GP flow stability analysis for N > 2 in (3) shall clarify the reason for such a discrepancy. 

Obtained results allow to fill a gap in the non-linear theory [7, 8] in which up to now instead 
of the linear exponential instability it was necessary to consider a stage of the seed algebraic 
instability (where small initial disturbance may grow only locally in time tending to zero with 

Noted herein new mechanism of the linear instability of the GP flow (taking into account 
possibility of existence of the bi-periodic initial disturbances described by the complementing to 
the Reynolds number Re parameter p being the ratio of two longitudinal periods) may be useful 
for the flow in a duct of square cross-section [26], the plane Poiseuille and Couette flows 
stability analysis. For the latter flows, as for the GP flow, also currently there is no matching of 
the experimental data and conclusions of the linear stability theory taking into account only pure 
periodic disturbances. For example, in relation with noted in [27] linear destabilization effect for 
the plane Poiseuille flow (when "the region where solid-fluid coexist stretches towards the wall" 
[27, p. 33]), pointed out bi-periodicity of the initial disturbance can, as in [12] (see footnote 2 on 
the page 1), be defined by the complementing character scale of the solid units transported by the 
flow of the fluid. 

We are grateful to S.I. Anisimov, E. A. Novikov and NA. Inogamov for useful comments 
and interest to the work. 

Appendix A 

1. From (3) and (5) for TV = 2 in Galerkin's approximation (representing disturbance 
dependency on the longitudinal coordinate z when generally a ^ /? in (5)), one gets for 

A = Aj + lA 2 '. 

V 2 

A i =- r 2 i -47r 2 fi 2 p 2 -^(a 1 ±j=D( ) , (Al) 

V 2 

A 2 = -2nj3pP n Re- i(a 2 ±\^ ) , (A.2) 

1/2 1/2 

/ / 2 2 2 

where D.=d' + 1 , D~=d' - 1 , I = a^ - a_ + 4^ Re 1 

a \ =Y \,2 ~ Y \,\ +47i 2 (3 2 {\-p 2 ), a 2 =27fiRe(P 22 ~pP n ), 

2 Ro 2 P P p 2 /3 2 2 2 2 

Re 2 = zi iz ^ ^ =l z +4(a l a 2 -2Re^d l ) z , ^=-45, 

(1-P) 



c. = -4Sctg7r(p H — ) , and S is defined in (6). 
1 P 

Condition X x > leads to the inequality 

(aRe+b) 2 >c + — T 
Re 2 

AP P P 2 p 2 S 
a = \\ 2 >b = nP(P 22 - P P n )a v 

(1-/0 

^ = «3(^ 1 +4^ 2 /?V)(^ 2 +4^ 2 /? 2 ) , « 3 =7 2 2 +/ 2 1 +4^ 2 /? 2 (l + ^ 2 ), 

??? 9 1Z Zl r, 

c = a 2 /? 2 (;r 2 (P -/?P ) 2 -4 P-). 

J (l-/?) 2 

2. In the limit of Re » 1 , inequality (A.3) for c > is reduced to the inequality 



where 



4^-b- 



Rc>Rfi tt tf) = . 



(-/ 



In (A.4), function Re , (/?) has the minimal value (given in (6)) when 



(A.3) 



(A.4) 



'->o-£ 



A(l + p 2 ) + B(l-p 2 ) 



where A and 5 are defined in the main text (see (6)) for A > 
3. On the neutral curve 

velocity V p I V m is as follows 



(A.5) 



3. On the neutral curve with 1 - (i.e. when equality takes place in (A.4)), the phase 



VJV m =pP u +P 22 ±(- 



LJ 2 



")/ 



(A.6) 



where (3 = /? 12 (Re, p) corresponds to the replacement of inequality by equality in (A.4). With 

such replacing of inequality by equality in (A.4), one gets quadratic equation with respect to /? . 
Its solution is as follows 



P = fl. 



Re±7Re 2 -Re 



In' S, 



(A.7) 



with Re > Re (/l , Re rt is from (6), when /? 12 =/? for 



Re = Re rt , and 



^ 1 =^^ 2 + i)Vr^4^-a-/)^L) l ^ l(1 2 p) 



Sb l {\S\p 2 P, 2 P 2 , 



P 2 P 2 A 2 



b x =P 22 - P P U ,S = - * '£" ctgrr(p + -). 

(1-/7) 71 b x 



1 
P 
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Fig. 1. Family of the curves of the neutral stability (with /L =0) according to (6) and (A.5). In 

Fig. la), a fragment from Fig. lb) is shown. In Fig. la), the upper curve (with two branches) 
corresponds to /? =0.463 (for p= 1.527), and the lower one to /?„ =1.099 (for/?=2.239) according 
to (A.5). Fig. la) is overlapped with a figure from [25] (see Fig. 12 in [25]) formally using l/2p= 

aS . In [25], a is the disturbance wave number, and 8 is the thickness of the boundary layer 
shift when streamlining a thin plate. In Fig. la), dots and dashed lines correspond to the 
experiment [25], and solid lines correspond to the theory of Schlichting (lower one with the 
Roman digit I), and Lin's (the upper one with the Roman digit II). 

5. Table of values of Re.^ and p obtained by the formulae (A.3) and (6) for p corresponding 
to the local minima of Re f , in the approximate formula (6) 



p 


0o (from 
(A.3)) 


Re, A (from 
(A.3)) 


0o (from 
(A.5)) 


Re th ( from 
(6)) 


1,527 


0,471 


448,455 


0,463 


442,278 


0,674 


1,124 


680,307 


1,101 


678,482 


0,447 


1,368 


1095,455 


1,358 


1093,824 


2,239 


1,100 


1983,171 


1,099 


1981,838 


2,791 


0,220 


13095,398 


0,219 


13095,285 


0,359 


1,114 


23816,499 


1,114 


23816,488 



J ■:-;- : 


2 s^~ 


'5 ' f~ " "'"" 


A 




-■■j 



Re 



0,807- 
0,807- 
0,806- 
0,806- 
0,805- 
0,805- 
0,804- 
0,804- 
0,803- 
0,803- 
0,802- 
0,802- 
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Fig. 2, a) Results of the experiments from [22] ] (blue triangles are used for the turbulent spot 
rear front velocity, and white triangles reflect velocity of the leading front) and [24] (squares and 
rounds for the velocity of the rear front). Phase velocity of the wave solutions [8] (numbers 1-5 
denote a degree of the azimuthal symmetry of the running wave). The velocities are scaled by the 
mean flow velocity in the pipe. Phase velocity (A.6) (upper "straight line" corresponds to the 
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plus sign in (A.6), and the lower one corresponds to minus) for p«l.53, p =0,471 that 
corresponds to the absolute minimum Re ,. =448.5 according to (A.3). b) Enlarged 

representation of the lower "straight line" from Fig. 2, a) according to (A.6) (the upper branch in 
Fig. 2, b) corresponds to the plus sign in (A.7), and the lower one to minus. 

Appendix B 

At first, let's consider conventional disturbance representation [2, 3]. Disturbance of the velocity 
field V(p , described by (1), may be represented as the following two independent modes: 

/L r + Htwcx N r 

V m = V , = e * I a JJy, -) 

, (B.l) 

/Lr + Unfix N r 

V <pl , n l v/ \,n r 

n = \ 

where a ,b are constants calculated by the use of the Galerkin's approximation for determining 

dependencies V , and V ~ on the radial coordinate r , x - — , z is the longitudinal coordinate 
along the axis of the pipe of radius R , J. is the Bessel function of the first order ( 
•A (V-\ ) — 0,m- 1,2,..). For such disturbance representation, it is known [2, 3], that real parts 
of the orders of exponential temporal evolution A and /L , defined independently for each mode 

out of these two, are negative, and disturbances tend to zero when r — > <x> for any Reynolds 
number. Obviously, it is true for the linear superposition of V x and V f2 . 

In our problem stating, for example, in the case of TV = 2 in (3), we consider already quite 
different from (B.l) representation (5) of the disturbance V in the form of the set of two now 

synchronous modes, i.e. the modes have one and the same order of exponential temporal 
evolution X (and also different periods along the longitudinal axis z for different radial modes): 

V<p = /'( V^O-u j) + A 20 e i2 ^ X Ur, 2 "-)) . (B.2) 

Here A may be defined only with the help of Galerkin's approximation for the used in (B.2) 
representation of the dependency of Vm on the longitudinal coordinate z , contrary to the exactly 

defined A and /L in (B.l). 

In Section 3 of the paper, we show that for the above-critical Reynolds number (6) real part of 
A in (B.2) can be already positive, i.e. GP flow can be exponentially unstable in the linear 
approximation. 
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